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Abstract. This paper is concerned with properties of maximal solutions of the 
Ricci and cross curvature flows on locally homogeneous three-manifolds of type 
SL2(IR). We prove that, gcncrically, a maximal solution originates at a sub-Riemannian 
geometry of Heisenberg type. This solves a problem left open in earlier work by 
two of the authors. 



1. Introduction 

1.1. Homogeneous evolution equations. On a closed 3-dimensional Riemannian 
manifold {M,g), let Rc be the Ricci tensor and R be the scalar curvature. The 
celebrated Ricci flow |Ham82j starting from a metric Qq is the solution of 



:i.ii 



i9 = -^Rc 

9(0) = 90- 



Another tensor, the cross curvature tensor, call it h, is used in |CH04] to define the 
cross curvature flow (XCF) on 3-manifolds with either negative sectional curvature 
or positive sectional curvature. In the case of negative sectional curvature, the flow 
(XCF) starting from a metric Qq is the solution of 



:i.2) 



= -^h 
9(0) = go- 



Assume that computations are done in an orthonormal frame where the Ricci tensor 
is diagonal. Then the cross curvature tensor is diagonal and if the principal sectional 
curvatures are fci,/c2;^3 {ki = Kj^j^, circularly) and the Ricci and cross curvature 
tensors are given by 

(1.3) Rii = kj + ki (circularly in I), 
and 

(1.4) ha = kjki. 
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A very special case arises when the 3-manifold (M, g) is locally homogeneous. In 
this case, both flows reduce to ODE systems. As a consequence, the cross curvature 
flow can be deflned even if the sectional curvatures do not have a deflnite sign, as this 
is the case for most of the homogeneous 3-manifolds. 

We say a Riemannian manifold (M, g) is locally homogenous if, for every two points 
X G M and y G M, there exist neighborhoods U of x and V of y, and an isometry 
(fi from {U,g\u) to {V,g\v) with ip{x) = y. We say that {M,g) is homogeneous if 
the isometry group acts transitively on M, i.e., U = V = M for all x and y. By 
a result of Singer |Sin60] , the universal cover of a locally homogeneous manifold is 
homogeneous. 

For a closed 3-dimensional Riemannian manifold (M, go) that is locally homoge- 
neous, there are 9 possibilities for the universal cover. They can be labelled by the 
minimal isometry group that acts transitively: 

(a) H{3) {H{n) denotes the isometry group of hyperbohc n-space); SO (3) x M; 
H{2) X M; 

(b) M?- SU(2); SL(2,M); Heisenberg; 1) = Sol (the group of isometry of plane 
with a flat Lorentz metric); £'(2) (the group of isometries of the Euclidean 
plane). This is called the Bianchi case in |IJ92j . 

The crucial difference between cases (a) and (b) above is that, in case (b), the universal 
cover of the corresponding closed 3-manifold is (essentially) the minimal transitive 
group of isometries itself (with the caveat that both SL(2,M) and E{2) should be 
replaced by their universal cover) whereas in case (a) this minimal group is of higher 
dimension. Case (b) corresponds exactly to the classiflcation of 3-dimensional simply 
connected unimodular Lie groups (non-unimodular Lie groups cannot cover a closed 
manifold). Any (M, g^) in case (b) is of form M = G/H, where G = M is the universal 
cover, if is a co-compact discrete subgroup of G, and the metric go descends from a 
left-invariant metric go on G. 

The cases of H{3), S0(3) x M, H{2) x M and all lead to well-understood and 
essentially trivial behaviors for both the Ricci and cross curvature flows. 

The forward behavior of the Ricci flow on locally homogeneous closed 3-manifolds 
was flrst analyzed by Isenberg and Jackson |IJ92j . The forward and backward be- 
haviors of the cross curvature flows are treated in |CNSC08| ICSCOSj whereas the 
backward behavior of the Ricci flow is studied in |CSC09j . Related works include 
[GliOSl IKMOll ILot07] . In |CSC09l ICSCOSj . the following interesting asymptotic be- 
havior of these Ricci and cross curvature flows in the backward direction was observed: 
Let gt be a maximal solution deflned on {—Tb,Tf) and passing through a generic go 
at t = 0. Then either Tb = oo and g{t) = e^^go, or < oo and there is a positive 
function r{t) such that r'^{t)g{t) converges to a sub- Riemannian metric of Heisen- 
berg type, see |Mon02l ICSCOSj . More precisely, in |CSC09I ICSCOSj . this result was 
proved for all locally homogeneous closed 3-manifolds, except those of type SL(2, M). 
Indeed, the structure of the corresponding ODE systems turns out to be somewhat 
more complicated in the SL(2,]R) case. 
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The aim of this paper is to prove the result described above in the case of locally 

homogeneous 3-manifolds of type SL(2,R), i.e., SL(2,M). This will finish the proof 
of the following statement announced in [CSC081 ICSCOQj . 

Theorem 1.1. Let {M,go) be a complete locally homogeneous 3-manifold (compact 
or not), corresponding to the case (h) discussed above. Let g{t),t G {—Th,Tf), be the 
maximal solution of either the normalized Ricci flow fl2.1l) or the cross curvature flow 
(11. 2p passing through go at t = 0. Let d(t) be the corresponding distance function on 
M. Assume that go is generic among all locally homogeneous metrics on M . Then 

• either M is of type , = oo and g{t) = go, 

• or Ti) < oo and there exists a function r{t) : {—Tf,, 0) — )■ (0, oo) such that, as 
t tends to —Tb, the metric spaces {M,r(t)d(t)) converge uniformly to a sub- 
Riemannian metric space (M, d{Tb)) whose tangent cone at any m & M is the 
Heisenberg group H3 equipped with its natural sub-Riemannian metric. 

Each of the manifolds M considered in the above Theorem is of the type M = 
G/H, where G is a simply connected unimodular 3-dimensional Lie group and H 
is a discrete subgroup of G. By a locally homogeneous metric on M, we mean a 
metric that descends from an invariant metric on G. A generic locally homogeneous 
metric on M is a metric that descends from a generic invariant metric on G. Note 
that homogeneous metrics on G can be smoothly parametrized by an open set f2 in a 
finite dimensional vector space and generic can be taken to mean "for an open dense 
subset of Q" . 

By definition, the uniform convergence of metric spaces (M, dt) to (M, d) means 
the uniform convergence over compact sets of {x,y) — )■ dt(x,y) to {x,y) — )■ d{x,y), 
see [CSCOSl ICSCOQj for more details. 

1.2. The Ricci and cross curvature flows on homogeneous 3-manifoIds. As- 
sume that is a 3-dimensional real Lie unimodular algebra equipped with an oriented 
Euclidean structure. According to J. Milnor |Mil76] there exists a (positively oriented) 
orthonormal basis (61,62,63) and reals Ai,A2,A3 such that the bracket operation of 
the Lie algebra has the form 

[ci, Cj] = Afc6fc (circularly in i,j, k). 

Milnor shows that such a basis diagonalizes the Ricci tensor and thus also the cross 
curvature tensor. If /j = ajakCi with nonzero ai,aj,ak G M, then [fi,fj] = XkCiifk 
(circularly in i,j, k). Using a choice of orientation, we may assume that at most one 
of the Xi is negative and then, the Lie algebra structure is entirely determined by 
the signs (in {— 1,0, +1}) of Ai,A2,A3. For instance, -|-, -|-, -|- corresponds to SU(2) 
whereas +, +, — corresponds to SL(2,]R). 

In each case, let e = (ei, 62, 63) G { — 1, 0, +1}^ be the corresponding choice of signs. 
Then, given e and an Euclidean metric go on the corresponding Lie algebra, we can 
choose a basis /i, /2, fs (with fi collinear to 6j above) such that 



(1.5) 



[fi: fj] = '^^kfk (circularly in i, j, k) . 
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We call a Milnor frame for qq. The metric, the Ricci tensor and the cross cur- 
vature tensor are diagonalized in this basis and this property is obviously maintained 
throughout either the Ricci flow or cross curvature flow. If we let be the dual 
frame of the metric go has the form 

(1.6) go = Aof ® + Bof ®f + Cof ® f. 

Assuming existence of the flow g{t) starting from go, under either the Ricci flow or 
the cross curvature flow (positive or negative), the original frame {fi)f stays a Milnor 
frame for g(t) along the flow and g{t) has the form 

(1.7) g{t) = A{t)f ®f + B{t)f ®f + C{t)f ® f. 

It follows that these flows reduce to ODEs in {A,B,C). Given a flow, the explicit 
form of the ODE depends on the underlying Lie algebra structure. With the help of 
the curvature computations done by Milnor in |Mil76] , one can find the explicit form 
of the equations for each Lie algebra structure. The Ricci flow case was treated in 
jIJ92] . The computations of the ODEs corresponding to the cross curvature flow are 
presented in [CNSC081 [CSC08] . 

1.3. Invariant metrics on SL(2,M). This paper is devoted to study of the Ricci 
and cross curvature flows on 3-dimensional Riemannian manifolds that are covered by 

SL(2,]R). Since it makes no differences, we focus on SL(2,M). Given a left-invariant 
metric go on SL(2, M), we fix a Milnor frame {fi}\ such that 

[/2,/3] = -2/l, [/3,/i]=2/2, [/l,/2]=2/3 

and 

go = Aof ®f + Bof ®f + Cof ® f. 
The sectional curvatures are 

K{f2 A /3) = J^i-^A' + B^ + C'- 2BC - 2AC - 2AB), 
K{h A /i) = ]4^(-35' + A^ + C^ + 2BC + 2AC - 2AB), 

KUi A /2) = t477(-3C" + ^2 + 52 + 2BC - 2 AC + 2AB). 
ABG 

Recall that the Lie algebra sl(2, M) of SL(2, M) can be realized as the space of two 
by two real matrices with trace 0. A basis of this space is 

These satisfy 

[H, V] = -2W, [W, H] = 2 V, [V, W] = 2H. 

This means that {W, V, H) can be taken as a concrete representation of the above 
Milnor basis {fi, f2, fs)- In particular, /i corresponds to rotation in SL(2,]R). Note 
further that exchanging /2, /s and replacing /i by — /i produce another Milnor basis. 
This explains the B, C symmetry of the formulas above. 
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1.4. Normalizations. Let g{t),t G /, be a maximal solution of 

n 8) I m9 = -2^ 

where v denotes either the Ricci tensor Rc or the cross curvature tensor h. By 
renormalization of g{t), we mean a family g{t),t G /, obtained by a change of scale 
in space and a change of time, that is 

g(t) = moit), t = 4>{t) 

where is chosen appropriately. The choices of are different for the two flows 
because of their different structures. For the Ricci flow, take 

= / ilj{s)ds. 
Jo 

In the case of the cross curvature flow, take 

t 



(j){t) = / ilj'{s)ds. 
Jo 

Now, set ip{t) = il){t). Then we have 

dt \dt J 

where v is either the Ricci or the cross curvature tensor of 'g. 

On compact 3-manifolds, it is customary to take ^hiip = |iJ, where v = ^ ^j^J^'^^ 
is the average of the trace of either the Ricci or the cross curvature tensor. In both 
cases, this choice implies that the volume of the metric 'g is constant. Obviously, 
studying any of the normalized versions is equivalent to studying the original flow. 
Notice that the finiteness of Tf, or Tj is not preserved under different normalization 
of flows. 



2. The Ricci Flow on SL(2,R) 

2.1. The ODE system. Mostly for historical reasons, we will consider the normal- 
ized Ricci flow 

(2.1) ^ = -2Rc+'^Rg, giO) = go, 

where go is a left-invariant metric on SL(2, M). Let g{t), t G {—Tb, Tf) he the maximal 
solution of the normalized Ricci flow through go. In a Milnor frame {fi}f for go, we 
write (see (11. 7p ) 

g = Af'®f'+Bf^f + Cf^f. 



f dA 


2 


'dt 


~3 


dB 
{ — 

dt 


2 
~3 


dC 


2 


^ H 


~3 
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Under (12. ip . ABC = AqBqCq is constant, and we set AqBqCq = 4. For this normalized 
Ricci flow, A, B, C satisfy the equations 

d A 9 

=r.[-A\2A + B + C) + A{B- CY' 
(2.2) \ ^ =^[-B^(2B + A-C) + B{A + Cf' 

Bf] 

2.2. Asymptotic results. Because of natural symmetries, we can assume without 
loss of generality that -Bo > Co- Then _B > C as long as a solution exists. Throughout 
this section, we assume that Bq> Cq. 

Theorem 2.1 (Ricci flow, forward direction, |IJ92] ). The forward time Tf satisfies 
Tf = oo. As t tends to oo, B — C tends to exponentially fast and 

Bit) ~ (2/3)t, C{t) ~ (2/3)t and A{t) ~ 

In the backward direction, the following was proved in [CSC09] . 

Theorem 2.2 (Ricci flow, backward direction). We have Tb G (0, oo), i.e., the max- 
imal backward existence time is finite. Moreover, 

(1) // there is a time t < such that A{t) > B{t) then, as t tends to —Tb, 

Ait) ^viit + Tb)-"\ Bit) ~ r/2(t + TbY'\ Cit) ~ r^.it + Tbf'' 

with rji = -\/6/4 and constants rji E (0, oo), i = 2,3. 

(2) // there is a time t < such that A < B — C then, as t tends to —Tb, 

A{t) ~ r/i(t + Tbf'\ B{t) ~ r/2(t + Tb)-''\ C{t) ~ 7^,{t + Tbf" 

with r]2 = "\/6/4, and constants rji G (0, oo), i = 1, 3. 

(3) If for all time t<0, B — C<A<B then, as t tends to —Tb, 

Ait) ~ Bit) ~ ^(t + Tb)^'/\ Cit) ~ |(t + Tb). 

As far as the normalized Ricci flow is concerned, the goal of his paper is to show 
that the third case in the theorem above can only occur when the initial condition 
(740,-80,^*0) belongs to a two dimensional hypersurface. In particular, it does not 
occur for a generic initial metric qq on SL(2,R). 

Theorem 2.3. Let Q = {(a, 6, c) E M.^ : a > 0,b > c > 0}. There is an open 
dense subset Qq of Q such that, for any maximal solution git), t G i—Tb,Tf), of the 
normalized Ricci flow with initial condition (A(0), -B(O), C(0) E Qo, ast tends to —Tb, 

(1) ezther Ait) ~ (v^/4)(t + T,)-^^, Bit) ~ 7],it + Tb)^l\ Cit) ~ i^^it + Tb)^l^ 

(2) or Ait) ~ r/i(t + Tb)'/\ Bit) ~ (v^/4)(t + Tb)-^'\ Cit) ~ i^.it + Tb)^'' 
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In fact, let Qi (resp. Q2) be the set of initial conditions such that case (1) (resp. 
case (2) ) occurs. Then there exists a smooth embedded hypersurface Sq G Q such that 
Q11Q2 o,re the two connected components of Q\Sq. Moreover, for initial condition 
on So, the behavior is given by case (3) 0/ Theorem 12.21 

In order to prove this result, it suffices to study case (3) of Theorem 12.21 This is 
done in the next section by reducing the system (12. 2p to a 2-dimensional system. 

Remark 2.1. The study below shows that, when the initial condition varies, all values 
larger than 1 of the ratio ri2/rj^ are attained in case (1). Similarly, as the initial 
condition varies, all positive values of the ratio rjs/rji are attained in case (2). 

2.3. The two-dimensional ODE system for the Ricci flow. For convenience, 
we introduce the backward normalized Ricci flow, for which the ODE is 



(2.3) 



-A\2A + B + C) + A{B- Cf] 
-B\2B + A-C) + B{A + Cf] 
-C\2C + A- B) + C{A + Bf] 



dA _ 2 

rf^ ~ ~ 3' 
dB _ 2 

~ ~ 3' 
dC_ 2 

By Theorem 12.21 the maximal forward solution of this system is defined on [0, Tb) 
with Tfo < 00. 

We start with the obvious observation that if {A, B,C) is a solution, then 

t ^ {XA{Xh), \B{\h), XC{Xh)) 

is also a solution. By Theorem 12.21 there are solutions with initial values in Q such 
that B/A tends to 00 and others such that B/A tends to 0. Let Qi be the set of 
initial values in Q such that B/A tends to and Q2 be the set of those for which 
B/A tends to 00. Let Sq be the complement of Qi U Q2 in Q- Again, by Theorem 
12. 2[ for initial solution in Sq, B/A tends to 1. 

The sets Qi, Q2, Sq must be homogeneous cones, i.e., are preserved under dilations. 
Hence, they are determined by their projectivization on the plane A = 1. So we set 
b = B/A, c = C/A and compute 



(2.4) 



db/dt = 2A%{l + b){b-c - 1) 
dc/dt = -2A'^c{l + c){b-c+l). 



This means that, up to a monotone time change, the stereographic projection of any 
flow line of (12. 3p on the plane A = 1 is a flow line of the planar ODE system 

db/dt = b{l + b){b-c-l) 
dc/dt = -c(l + c)(6-c+ 1). 



(2.5) 



Set Q = {b > c > 0}. By Theorem 12. 21 any integral curve of (12. 5p tends in the forward 
time direction to either (0, 0), (00, Cqo), < Cqo < C)0, or (1, 0). The equilibrium points 
of (12. 5p . i.e., the points where {db/dt, dc/dt) = (0, 0), are (1, 0), (0, 0) (notice that they 
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are in dQ but not in Q). To investigate the nature of these equihbrium points, we 
compute the Jacobian of the right-hand side of the fl2.5p which is 

/ - 36c - c - 1 -b^-b \ 

\ -c^-c 3c2 - 26c - 6 - 1 J ■ 

In particular, at (0, 0), this is —1 times the identity matrix and the equihbrium point 

/ 2 -2 \ 

(0,0) is attractive. At (1,0), the Jacobian is ( g 2 ]' '^'^^^ point is a hyperbohc 

saddle point (the eigenvalues are 2 and —2). By Theorem 12.21 any integral curve of 
(12. 5 p ending at (1, 0) must stay in the region {6 — c < 1 < 6}. In that region, 6 and c 
are decreasing functions of time t and ^ is positive. Using this observation and the 
stable manifold theorem, we obtain a smooth increasing function : [1, oo) — )■ [0, oo) 
whose graph 7 = {(6, c) : c = 0(6)} is the stable manifold at (1, 0) in Vt. By Theorem 
12.11 and (12. 5p . 7 is asymptotic to c = 6 at infinity, and 0'(1) = 2. In particular, f2 \ 7 
has two components VLi and where (0,0) G i^i. Further any initial condition in 
VL whose integral curve tends to (1,0) must be on 7. It is now clear that the cases 
(1), (2) and (3) in Theorem 12.21 correspond respectively to initial conditions in fii, 
^2 and 7. This proves Theorem 12.31 with Qi the positive cone with base VLi and 5*0 
the positive cone with base 7. 

Figure ([T]) shows the curve 7 and some flow lines of (12. 5p . It is easy to see from 
(12. 5p that the flow lines to the right of 7 have horizontal asymptotes {c = Cqo} and 
that all positive values of Coo appear. This proves Remark 12.1 1 in case (2) of Theorem 
12.31 The proof in case (1) is similar, but a different choice of coordinates must be 
made. 




1 2 3 4 5 6 



Figure 1. The flow line diagram of system (12. 5p 
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3. The cross curvature flow on SL(2,M) 

3.1. The ODE system. We now consider the cross curvature flow (11. 2p where qq is 
a left-invariant metric on SL(2,M). Let g(t), t G {—Tf,,Tf), be the maximal solution 
of the cross curvature flow through Qq. Writing 

g = Af®f' + Bf®f + Cf^f, 

we obtain the system 



(3.1) 



' dA 


2AF2FS 


H ~ 


{ABCy 


dB 




'dt ~ 


(ABC)^ 


dC 


2CF1F2 


^ Itt ~ 


{ABCf 



where 

Fi = - 3^2 + ^2 + - 2BC - 2AC - 2AB, 
F2 = - 3B^ + A'^ + C^ + 2BC + 2AC - 2AB, 
F3 = - 3C^ + A'^ + B'^ + 2BC - 2AC + 2AB. 

3.2. Asymptotic Results. Without loss of generality we assume throughout this 
section that Bq > Co- Then B > C as long as a solution exists. The behavior of the 
flow in the forward direction can be summarized as follows. See |CNSC08] . 

Theorem 3.1 ((XCF) Forward direction). If Bq = Co then Tj = 00 and there exists 
a constant A^o G (0, 00) such that, 

B{t) = C{t) ~ (24Aoot)^/^ and A{t) ^ A^ as t ^ 00. 

If Bq > Co, Tf is finite and there exists a constant E G (0, 00) such that, 

A{t) ~ Bit) ~ E{Tf - ty^/^ and C{t) ~ 8(Tj - tf'^ as t Tf. 

In the backward direction, the following was proved in [CSCOSj . 

Theorem 3.2 ((XCF) Backward direction). We have Tf, G (0, 00), i.e., the maximal 
backward existence time is finite. Moreover, 

(1) // there is a time t < such that A(t) > B(t) — C(t) then, as t tends to —Th, 

A{t) ^viit + nr'^'\ Bit) ~ ri^it + ^,)3/l^ Cit) ~ r,,{t + T,)^/'' 

for some constants rji G (0, 00), z = 1, 2, 3. 

(2) // there is a time t < such that A < \ {2VB^ - BC + - B - C) then, 
as t tends to —Tb, 

A{t) ^viit + n)'/'\ Bit) ~ r^2it + ^,)-l/l^ Cit) ~ r/3(t + n)'/'' 

for some constants rji G (0, 00), i = 1, 2, 3. 
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(3) If for all time t < 0, | (2^52 - BC + -B-C)<A<B-C then, as t 
tends to —Th, 

A{t) ~ QAr]~\t + Tf,), B{t) ~ 77 + A^y2{t + n), C{t) ^ - A^2{t + T,) 
for some rj G (0, 00). 

As far as the cross curvature flow is concerned, tlie goal of this paper is to show 
that the third case in the theorem above can only occur when the initial condition 
{Aq, Bq,Cq) belongs to a two dimensional hypersurface. In particular, it does not 
occur for a generic initial metric go on SL(2,R). 

Theorem 3.3. Let Q = {{a,b,c) G : a > 0, 6 > c> 0}. There is an open dense 
subset Qq of Q such that, for any maximal solution g{t), t G (— Tf,,Tj), of the cross 
curvature flow with initial condition {A{0), -B(O), C(0)) G Qo, as t tends to —T^, 

(1) either A{t) ^ Vi{t + B{t) ~ r^^it + T.f/^^ C{t) ~ ^t + T.f/^^ 

(2) or A{t) ~ r/i(t + T,)3/l^ Bit) ~ //^(t + T,)-Vi^ ~ Vsit + 

In fact, let Qi (resp. Q2) he the set of initial conditions such case (1) (resp. case [2)) 
occurs. Then there exists a smooth embedded hypersurface So C Q such that Qi,Q2 
are the two connected components of Q\So- Moreover, for initial condition on So, 
the behavior is given by case (3) of Theorem 13.21 

In order to prove this result, it suffices to study case (3) of Theorem 13.21 In that 
case, it is proved in |CSC08] that A, B and C are monotone {A, B non- decreasing, C 
non-increasing) on {—Tb, 0]. In order to understand the behavior of the solution, and 
because of the homogeneous structure of the ODE system (13. ip . we can pass to the 
affine coordinates {A/C, B /C). This leads to a two-dimensional ODE system whose 
orbit structure can be analyzed. 

Remark 3.1. The analysis below shows that, in case (1) of Theorem \3.3\ and when 
the initial condition varies, all the values larger than 1 of the ratio ri2/ri^ are attained. 
Similarly, in case (2), all the values of the ratio rji/rj^ are attained. 

3.3. The two-dimensional ODE system for the cross curvature flow. For 

convenience, we introduce the backward cross curvature flow, for which the ODE is 

' dA _ 2AF2F3 
H ~{ABCf' 
\ dB _ 2BF3F1 

^^■^^ I 'dt "(ABCr ' 

dC _ 2CF1F2 
^ lit ~{ABCy' 

where {Fi}l are deflned as before. By Theorem 13.21 the maximal forward solution of 
this system is deflned on [0,Tfe) with < 00. 
Note that if {A, B, C) is a solution, then 

t ^ {XA{t/X^), XB{t/X^), XC{t/X^)) 
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is also a solution. By Theorem 13. 2[ there are solutions with initial values in Q such 
that A/B tends to oo and others such that {A/B,C/B) tends to (0,0). Let Qi be 
the set of initial values in Q such that A/B tends to oo and Q2 be the set of those 
for which {A/B,C/B) tends to (0,0). Let 5*0 be the complement of Qi U Q2 in Q. 
Again, by Theorem 13.21 for initial solution in So, {A/B,C/B) tends to (0, 1). 

The sets Qi, Q2, So must be homogeneous cones, i.e., are preserved under dilations. 
Hence, they are determined by their projectivization on the plane B = 1. So we set 
a = A/B, c = C/B and compute 



(3.3) 
where 




= + l)(« + c- 1)03 

= -(Bl^c(l-c)(a + c+l)0i, 

3a^ + 1 + - 2c - 2ac - 2a, 
3c^ + + 1 + 2c - 2ac + 2a. 



(3.4) 



- 

This means that, up to a monotone time change, the stereographic projection of 
any flow line of (13. 2p on the plane i? = 1 is a flow line of the planar ODE system 

da/dt = a(a + l)(a + c — 1)03 
dc/dt = — c(l — c)(a + c + 1)01. 

Set r2 = {a>0,l>c>0}. By Theorem 13. 2[ any integral curve of (13. 3p tends 
in the forward time direction to either (0,0), (oo,Coo), < Coo < 00, or (0,1). 
The equilibrium points of (13. 3p . i.e., the points where {da/dt, dc/dt) = (0,0), are 
(0, 0), (1, 0), (0, 1). To investigate the nature of these equilibrium points, we compute 
the Jacobian of the right-hand side of the ODE which is 

where 

Yn = (3a^ + 2ac + c- 1)03 + 2a(a + l)(a + c- l)(a - c+ 1), 

Y12 = a(a + l)[03 + 2(a + c-l)(-3c+l-a)], 

Y21 = -c(l -c)[0i -2(a + c+l)(3a + c+l)], 

Y22 = (3c2 + 2ac-a-l)0i-2c(l-c)(a + c+l)(c-l-a). 

In particular, at (0,0), this is —1 times the identity matrix and the equilibrium 
point (0, 0) is attractive. Its basin of attraction corresponds to region Q2 defined 

above. At (1,0), the Jacobian is ^ g 8 ^' "^^^^ point is a repelling fixed point. It 

reflects the behavior of the forward Ricci flow described in Theorem \3.1i B > C). At 

(0, 1), the Jacobian is ( n ) • This is the equilibrium point of interest to us and a 



^00^ 

more detailed analysis is required to determine trajectories that tend toward it. This 
is done with coordinate transformations that blow-up the equilibrium. 
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(3.5) 



(3.6) 



Blow-up transformations introduce coordinates in which the blown up equilibrium 
becomes a circle or projective line representing directions through the equilibrium. 
Blow-up transformations reduce the analysis of flows near degenerate equilibria to 
flows with less degenerate equilibria. The blown up system allows us to analyze 
the trajectories that are asymptotic to the equilibrium. In particular, trajectories 
approaching the equilibrium of the original system from different directions yield 
different equilibria in the blown up system. 

Translating the equilibrium point to the origin by setting e = c — 1, the equations 
(13. 3p become 

da/dt = -a(a + l)(a + e)(3e2 + 4e + 2ae - a^), 
de/dt = e(e + l)(a + e + 2)(e^-2ae-4a-3a2). 

The leading order terms of equations (13.51) have degrees 3 and 2: the next coordinate 
transformation a = u'^,e = v oi the region a > produces a system in which the 
leading terms of both equations have degree 3: 

du/dt = -lu{u'^ + l){v + u^){3v'^ + Av + 2u^v-u'^), 
dv/dt = v{v + l){v + 2 + u'^){v'^ -2u'^v-4u'^ -Su"^). 

To blow up the origin, the equations (13. 6 p are transformed to polar coordinates 
{u,v) = (r cos(^), r sin(^)) and then rescaled by a common factor of r, yielding the 
vector field X defined by 
(3.7) 

dr/dt = l/2r[-33r3(cos(^))^(sin(^))^ -29r2(cos(^))^(sin(^))^ 

-11 (cos {e)f{sm{e)f — (cos (^))^sin (6) - 5 (cos (^))'sin {9) 

(cos {9)Y^ - 35 r (cos {9))^ (sin {9)f - 20 (cos {9)f (sin {9)f 
— 4r cos (9))^ sin {9) + r^ {cos{9)f + 2 (sin {9)) 
-2 r3 (sin {9)f (cos {9)f - 18 (cos {9)f (sin (9))^ 
-10 (sin (e))^ (cos (9))^ + 6 r (sin {9)f 
-6 (sin {9)f (cos {9)f + 4 (sin (^))^], 
d9/dt = -l/2cos{9)sm{9)[-2r^sm{9))'^-r^sm{9)f{cos{9)f 

+9r2 (cos (^))^(sin {9)f + 5 (sin (9))'^ (cos (9))'^ — 9 r (sin (9))^ 
-f28r(cos(e))^sin (9) + 25 (cos (6'))^ sin (9) + 5 sin (9) (cos (6'))^ 
-8 (sin {9)f + 16 (cos {9)f + 20 (cos (9))^ 
+7r^ (cos {9)f + r^ {cos{9)f]. 

In these equations, the origin of equations (13. 6 p is blown up to the invariant circle 
r = 0, and the complement of the origin becomes the cylinder r > 0. Trajectories 
that tend to the origin in equation (13. 6p yield trajectories that tend to an equilibrium 
point of equations (13. 7p on the circle r = 0. Now the zeros of d9/dt on the circle 
r = are equilibria of the rescaled equations obtained from (13. 7p . They are located 
at points where cos^(0) = 0, 1 or 1/3. The equilibria determine the directions in which 
trajectories of (13. 6 p can approach or leave the origin. These directions correspond to 
different approach to the equilibrium (0,1) of (13. 4p . In the (a, c) coordinates, the 
directions 9 = ±7r/2 correspond to approaching (0,1) along curves tangent to the 
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c-axis with tangency degree greater than 2. The direction ^ = corresponds to 
approaching (0, 1) along curves tangent to c = 1. The directions ±6q with cos^6'o = 
1/3 correspond to approaching (0,1) along curves asymptotic to the parabola a = 
|(c — 1)^. Observe that this is consistent with Theorem 13.2( 3). Our goal is to show 
that this can only happen along a particular curve. 

Since the circle r = is invariant, the Jacobians at the equilibria discussed above are 
triangular. The stability of each equilibrium is determined by the signs of ^{dr/dt) 
and d{d6/dt)/d6 when these are non-zero. The equilibria with cos^(6') = have 
^{dr/dt) = 2 and d{d9/dt)/d9 = —4, so the point is a saddle with an unstable 
manifold in the region r > 0. Equilibria with cos^(6') = 1/3 have ^{dr/dt) = —4/3 
and d{d9 / dt) / 89 = 16/3, so these equilibria are also saddles but with stable manifolds 
in the region r > 0. After change of coordinates, only one of these stable manifolds, 
call it 7o, belongs to the region {c < 1}. This curve 70 provides the only way to 
approach (0, 1) which is consistent with case (3) of Theorem 13.21 The equilibria with 
cos^(6') = 1 have ^{dr/dt) = and d{d9 / dt) / 39 = —8, so further analysis is required 
to determine the properties of nearby trajectories. However, because of Theorem 13.21 
it is clear that any solution of (13. 4p approaching the line {c = 1} has a — )■ 00, hence 
cannot approach (1, 0). The following argument recovers this fact directly from (13. 7p . 

Note that when 9 = 0, d9/dt = and dr/dt = (r^ + r'^)/2. Therefore, the r axis 
is invariant and weakly unstable. The trajectory along this axis approaches (0, 0) as 
t — )■ —00. To prove that no other trajectories in r > approach the origin as t — ±00, 
we consider the vector field Y defined by subtracting cos^(6')/2 + r^cos'^°(6')/2 from 
dr/dt in X: 



(3.8) 

'dr/dt = 1/2 sin (9) r[-5r^ {cos (9) f -4r {cos (9))^ 

-29 r2 (cos {9))^ sin {9) - 20 sin {9) (cos {9)f - 35 r (cos {9)f (sin {9)f 
-r^ (cos {9)Y - 11 (cos {9)f sin {9) - 33 (cos {9))'^ (sin {9)f 
-6 (cos {9)f (sin {9)f - 10 (cos {9))^ (sin {9)f 
-2 r3 (cos {9) f (sin (^))^ - 18 (cos {9) f (sin {9))^ + 6 r (sin (^))^ 
+2r2(sin(^))^ + 4 (sin(^))% 
d9/dt = -1/2 cos{9)sm{9)[-2r'^{sm{9))'^ -r^sm{9)f{cos{9)f 

+9 (cos {9)f (sin {9)f + 5 (sin {9)f (cos {9))^ - 9 r (sin {9)f 
+28 r (cos {9)f sin {9) + 25 (cos {9))'^ sin {9) 

+5 sin {9) (cos {9)f - 8 (sin {9)f + 16 (cos {9)f + 20 (cos {9))^ 
+7r4 (cos(^))' + r6 (cos {9))% 



The vector field Y is transverse to the vector field X in the interior of the first 
quadrant: d9/dt < for both X and Y and the r component of X is larger than 
the r component of Y. Therefore, trajectories of X cross the trajectories of Y from 
below to above as they move left in the {9, r) plane. The vector field Y has a common 
factor of sin(^) in its two equations. When Y is rescaled by dividing by this factor, 
the result is a vector field that does not vanish in a neighborhood of the origin. Since 
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the 9 axis is invariant for the Y trajectory 7 starting at (6',r), r > approaches 
the r axis at a point with r > 0. The X trajectory starting at (6*, r) hes above 7, so 
it does not approach the origin. This proves that the only trajectories of equations 
asymptotic to the origin he on the r and 9 axes. 
In conclusion, the above analysis shows that the regions Qi, Q2 are separated by 
the 2- dimensional cone 5*0 determined by the curve 70 in the (a, c)-plane. This proves 
Theorem 13.31 Figure |2] describes the flow lines of (13. 4 p . The part of interest to us 
is the part below the line {c = 1}, which corresponds to {B > C}. The part above 
the line {c = 1} corresponds to the case {B < C}, where the role of B and C are 
exchanged. The most important component of this diagram are the flow lines that 
are forward asymptotic to (0, 1). They correspond to the hypersurface Sq in Theorem 
f l3.3p . The flow lines in the upper-left corner have vertical asymptotes {a = a^} 
with all positive values of appearing. Similarly, the flow lines on the right have 
horizontal asymptotes {c = Coo} with all positive values of Coo appearing. These 
facts can easily be derived from the system (13. 4p . This proves the part of Remark 
13.11 dealing with case (1) of Theorem 13.31 The other case is similar using different 
coordinates. 




1 2 3 4 5 6 



Figure 2. The flow line diagram of system (13. 4p 
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